8 Coherent State Path Integral Quantization of 
Quantum Field Theory 

8.1 Coherent states and path integral quantization. 
8.1.1 Coherent States 

Let us consider a Hilbert space spanned by a complete set of harmonic oscillator 

states with n = 0, . . . , oo. Let and a bo a pair of creation and annihi- 

lation operators acting on that Hilbert space, and satisfying the commutation 
relations 

[a, at] = 1 , [at, at] = , [a, a] = (1) 
These operators generate the harmonic oscillators states {|n)} in the usual way. 

In) = -1= (at)" |0) (2) 



a\0) = (3) 



whore |0) is the vacuum state of the oscillator. 
Let us denote by \z) the coherent state 

\z) = e^«^ |0) (4) 
(^1 = (0| e^^ (5) 

where z is an arbitrary complex number and z is the complex conjugate. The 
coherent state \z) has the defining property of being a wave packet with opti- 
mal spread, i.e., the Heisenberg uncertainty inequality is an equality for these 
coherent states. 

How does a act on the coherent state \z)l 

oo „ 
Z 



Since 
we get 

Thus, we find 



a,(at)"] =n (at)"-' (7) 
«N) = E^ ([a, (at)"] + (at)" a) |0) (8) 



n=0 



7 —1 



n=0 



Therefore \z) is a right eigenvector of a and z is the (right) eigenvalue. 
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Likewise we get 



«'E^(«T 10) 

n=0 

oo „ 

oo 

n=0 

oo 



n — 



(n + 1)! 
(at)" |0) 



n=l 



Thus, 



(9z 



Another quantity of interest is the overlap of two coherent states, {z\z'), 

{z\z') = (0|e^^ e^'°^ |0) 

We will calculate this matrix element using the Baker-Hausdorff formulas 

1 



e-" = e 2 



A,B 



= e 



AS 



(10) 
(11) 

(12) 
(13) 



which holds provided the commutator 



A,B 



is a c- number, i.e., it is propor- 



tional to the identity operator. Since [d, at] — 1, we find 



But 
and 

Hence we get 



{z\z') = e^^' (0|e^'^^ e^^lO) 



e^'^ |0) = |0) 



(0| 



„z a 



(0| 



(14) 
(15) 
(16) 

(17) 



An arbitrary state lip) of this Hilbert space can be expanded in the harmonic 
oscillator basis states {|n)}, 



n=0 



E I-) = E !S i-r 10) 



n=0 



(18) 
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The projection of the state \ip) onto the coherent state \z) is 

{zW = f:^{z\{aT\0) (19) 

n=0 

Since 

(^1 a'' =z {z\ (20) 

we find 

n=0 

Therefore the projection of onto \z) is the anti-holomorphic {i.e., anti- 
analytic) function V'(^)- In other words, in this representation, the space of 
states {\ip)} are in one-to-one correspondence with the space of anti-analytic 
functions. 

In summary, the coherent states {l^)} satisfy 

a\z) — z\z) {z\a = dz{z\ 
a^z) = d^\z) {z\a^=z{z\ (22) 
(^|V)=V(^) m=ij{z) 



Next we will prove the resolution of identity 



dzdz 
2m 

Let IV') and be two arbitrary states 



/ 



e 



-ZZ\ 



z){z\ (23) 



OO , 

oo , 
n— 

(^1^) = E ^ (24) 

n=0 



Let us compute the matrix element 

(<^|/»=^^(n|/|m) (25) 

m.n 

Thus we need to find 

{n\i\m) = [ ^ e-\^\^{n\z){z\m) (26) 

J 2lTl 



Recall that the integration measure is defined to be given by 

dzdz dRezdlmz 
27ri TT 



(27) 
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where 

1 z" 
(n|^) = ^(0|(an.) = ^(0|z) (28) 
vn! vn! 

and 

{z\m) = -={z\ (at)™ |0) = ^{z\Q) (29) 
Now, since |(0|2;)p = 1, we get 

i^n\i\m)= [d_^e_^,n,r.^ r^,^ ^^^„+.,.(n - 
J 2m y/n\m\ Jo Jo 27r \/n\m\ 

(30) 

Thus, 

(n|/|m) = ^ / dx x"e~^ = (nlm) (31) 

TV- Jo 

Hence, we have found that 

= m) (32) 

for any pair of states \tp) and Therefore I is the identity operator in that 
space. We conclude that the set of coherent states {| 2:)} is an over-complete set 
of states. 

Furthermore, since 

(^1 (at)" (a)™ \z') = z''z'-^{z\z') = ^z'^^e^^' (33) 

we conclude that the matrix elements of any arbitrary normal ordered operator 
of the form 

i = ^^n,m(a^)"(a)'" (34) 

are equal to 

i,z\A\z') = (^A„,„^"2"" J e^^' (35) 

\n,m / 

Therefore, if A{d, a^) is an arbitrary normal ordered operator (relative to the 
state |0)), its matrix elements are given by 

(z|i(a,at)|z') ==^(z,z')e^^' (36) 

where A{z, z') is a function of two complex variables z and z' , obtained from A 
by the formal replacement 

a ^ z' , ^ z (37) 

For example, the matrix elements of the the operator N = a)a, which measures 
the number of excitations, is 

{z\N\z') = (^latal^') = zz' e^^' (38) 
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8.1.2 Path Integrals and Coherent States 

As usual we will want to compute the matrix elements of the evolution operator 
U, 

T - 
—i—H(aK a) 
U = e h ^ ' (39) 

where H{a^,a) is a normal ordered operator. Thus, if \i) and |/) denote two 
arbitrary initial and final states, we can write the matrix element of U as 



(/|e h ' ' ' 



However now, instead of inserting a complete set of states at each intermediate 
time tj (with j = 1,. . . ,N), we will insert an over-complete set at each 

time tj through the insertion of the resolution of the identity, 



{f\(l-i^H{a\a)y\i) 



-j-j- dzjdzj 



N 

■El 



=1 {f\(l-i^H{a\a))\z^) 
x{zn\ (l-i^H{a^,aj^ |^;jv-i) • • • (l - i^H {a^ , a)^ \i) 

N 

-El- 



dzjdzj 
2m 



N-l 



W {zk+i\ (l-i^H{d\a)) \zk) 



x(/| (l-^tH{d\d)) \zn){zA (l-i'-H{d\d)) \zi) 

In the limit e ^ these matrix elements are 

{zk+i\ {l-i^H{a\a)^ \zk) ={zk+i\zk) - i^{zk+i\H{a\d)\zk) 

={zk+i\zk) - ijH{zk+i,Zk] 



(41) 



(42) 



where H{zk+i, Zk) is a function which is obtained from the normal ordered 
Hamiltonian by the substitutions Zk+i and d ^ Zk- Hence, we can write 
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the following expression for the matrix element 

T - 

(/|e h '\i) = 



N 



N-1 



= lim 

e^O,Ar->oo 



/(n 

x{f\zN){zi\i) 



dzjdzj 
2m 



e eJ=i [] [ 



1 - i-H{zk+i,Zk) 



^ {zi\i) 



(43) 



By further expanding the initial and final states in coherent states 

dZfdZf 



(/I 



2iTi 



^fizf){zf\ 



dzidzi _|7.|2 , , , 



(44) 



we find 



T - 
—i—H{d\ a) 

{f\e ^ 'N) 



VzVz e 



dt 



2i 



{zdtz — zdtz) — H{z, z) 



(45) 



This is the coherent-state form of the path integral. We can identify in this 
expression the Lagrangian L as the quantity 



i = — {zdtZ - zdtz) - H{z, z) 



(46) 



Notice that the Lagrangian in the coherent-state representation is first order in 
time derivatives, because of this feature we are not guaranteed that the paths 
are necessarily differentiable. This property leads to all kinds of subtleties that 
for the most part we will ignore in what follows. 



8.1.3 Path integral for a gas of non-relativistic bosons at finite tem- 
perature. 

The field theoretic description of a gas of (spinless) non-relativistic bosons is 
given in terms of the creation and annihilation field operators (f>^{x) and 4>{x), 
which satisfy the equal time commutation relations (in d space dimensions) 

= S\x-y) (47) 
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Relative to the empty state |0), i.e., 

${x)\0) = (48) 
the normal ordered Hamiltonian is 



H= j (X^xPix) 

+\ j d^x j dS4>\x)^\y)u{s-y)ky)k^) 



(49) 



where m is the mass of the bosons, /i is the chemical potential, V{x) is an 
external potential and U{x — y) is the interaction potential between pairs of 
bosons. 

Following onr discussion of the coherent state path integral we see that it 
is immediate to write down a path integral for a thermodynamic system of 
bosons. The boson coherent states are now labelled by a complex Geld (j){x) and 
its complex conjugate ^{x). 



dx 6(x)<b{x) 

mx)}) = eJ |0) (50) 

which has the coherent state property 

^{x)m)=mm) (si) 

as well as the resolution of the identity 



r -[dx\cl>{x)f 
1= jv^Vct>*e J m)m\ (52) 

The matrix element between an initial state \i) and a final state |/) separated 
by a time span T = tf — U (not to be confused with the temperature!) now 
takes the form 

--HT 
{f\e h \i) = 

J V(j)V4> expS^^Jj' dt (^j d^'x y [(j){x, t)di4>{x, t) - 4>{x, t)dt(t>{x, t)] - H[4>, 

- ^ ,^ \ fdxm^,tf)\'' + \ct>{x,t,)\') 

x^f{<Pix,tf))^i{(j3{x,ti)) e^J 



where is 

i,2 



(53) 



H[,t>A=j d^xm -|;^v'+M + m m + ^Jd^x j d''y\ct>{x)?\4>i.y)?U{x-y) 

(54) 
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It is also possible to write the action S in the less symmetric but simpler form 

{dx = dtdx) 



j dx 



x) ( im + —V - At - V{x) ] ^{x) 



-\ j d'x J d''y\,^{x)\'\cp{y)\'Uix^y) 

(55) 

where U{x — y) = U{x — y)S{tx — ty). 

This formulation is useful to study superfluid Helium and similar problems. 
Suppose for instance that we want to compute the partition function Z for this 
system of bosons at finite temperature T, 

Z = tre~^^ (56) 

where /3 = 1/T (in units where ku = !)• The coherent-state path integral 
representation of the partition function is obtained by 

1. restricting the initial and final states to be the same \i) = \f) and arbitrary 

2. Summing over all possible states 

3. and finally a Wick rotation to imaginary time t —ir, with the time-span 
T — > —i/3h {i.e., Periodic Boundary Conditions in imaginary time) 

The result is the (imaginary time) path integral 

Z = J V(I)V4> e'^^^'^' ^) (57) 
where Se is the Euclidean action 

Se{4>,4>) = - j dr j dx ' 



-hdr+H - i^V^ -V{x) 
2m 



1 



\t I dx I dyU{x-y) \(t^{x)\'\^{y)\' 



(58) 



The fields (f){x) = (j){x, r) satisfy Periodic Boundary Conditions (PBC's) is imag- 
inary time 

^(f , r) = (j){x, T + ph) (59) 

This requirement suggests an expansion of the field ^{x) in Fourier modes of 
the form 

oo 

mr)= E e'^^^^m^^n) (60) 
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where the frequencies cj„ (the Matsubara frequencies) must be chosen so that 
(/) obeys the required PBCs. Wc find 

27^ 2ttT ^ , , 

w„ = -0a = -^'^ ' n e Z (61) 

where n is an arbitrary integer. 



8.2 Fermion Coherent States 

In this section we will develop a formalism for fermions which follows closely 
what we have done for bosons while accounting for the anti-commuting nature 
of fermionic operators, i. e. the Pauli Principle. 

Let {c|} be a set of fermion creation operators, with i = 1, . . . , TV, and {cj} 
the set of their N adjoint operators, i. e. the associated annihilation operators. 
The number operator for the i-th fermion is rii = c\ci. Let us define the kets 
|0i) and |lj), which obey the obvious definitions: 

c.|0.)=0 410.) = |1.) 

cU\0^)=Q 4c,|l,) = |l,) ^ ' 

For N fermions the Hilbert space is spanned by the anti-symmetrized states 
|ni, . . . ,njv). Let 

|0) = |0i,...,0jv) (63) 

and 

|ni,...,nAr) =cl...4|0) (64) 
As we saw before, the wave function (ni, . . . , tiatI^') is a Slater determinant. 



8.2.1 Definition of Fermion Coherent States 

We now define fermion coherent states. Let with i = 1, . . . , iV, be a 

set of 2N Grassmann variables (also known as the generators of a Grassmann 
algebra.) These variables satisfy, by definition, the following properties 

{6,0} = {CU,} = = = e7 = (65) 

We will also require that the Grassmann variables anti-commute with the fermion 
operators: 

{^i, Cj} = {^i, c\ } = {iu c„- } = [ii, ct } = (66) 
Let us define the fermion coherent states 

\0 ^ e-^^[\Q) (67) 
{i\ ^ (0|e^'^ (68) 
As a consequence of these definitions we have: 

e"^^^ = 1 - (69) 
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Similarly, if i/' is a Grassmann variable, we have 

(^IV;) = (Ojef^e-^'^^ |0) = 1 + = e^V' (70) 
For N fermions we have, 

JV 

|O = |a,---,M=n£ie-^^^l|0)^e -1 ^|0) (71) 
since the following commutator vanishes, 

6cU,c]]=0 (72) 
8.2.2 Analytic functions of Grassmann variables 

We will define V(0 be an analytic function of the Grassmann variable if it 
has a power series expansion in ^, 

V'(o = V'o + V'i^+V'2<e' + --- (73) 

where V'n € C. Since 

C = 0, Vn > (74) 

then, all analytic functions of a Grassmann variable reduce to a first degree 
polynomial, 

V(0 = V-o + V'lC (75) 
Similarly, we define complex conjugation by 



^(0 = Vo + V'l? (76) 

where -ipo and -ipi are the complex conjugates of tpo and ipi respectively. 
We can also define functions of two Grassmann variables $, and ^, 

= ao + aiC + aif + aiaf^ (77) 

where fli , 0.1 and ai2 are complex numbers; ai and ai are not necessarily complex 

conjugates of each other. 

1. Differentiation of Grassmann variables: 

Since analytic functions of Grassmann variables have such a simple struc- 
ture, differentiation is just as simple. Indeed, we define the derivative as 
the coefficient of the linear term 

= V-i (78) 

Likewise we also have 

d,m = ^i (79) 
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Clearly, using this rule we can write 

d^m = -dm) = (80) 

A similar argument shows that 

d^A{^,0 = ai-a,2^ (81) 

= ai+ai2^ (82) 

d^d^A{^,0 = -ai2 = -d^d^A{^,0 (83) 
from where we conclude that (9^ and anti- commute, 

{a^, 5j } = 0, and d^d^ = = (84) 

2. Integration over Grassm,ann variables 

The basic differentiation rule of Eq. (78) implies that 

1 = a^c (85) 

which suggests the following definitions: 



I 



I 



d^l = (86) 
d^d^^ = ("exact differential") (87) 
= 1 (88) 



/ 



Analogous rules also apply for the conjugate variables ^. 

It is instructive to compare the differentiation and integration rules: 

Thus, for Grassmann variables differentiation and integration are exactly 
equivalent 

d^<^ j dC (90) 
These rules imply that the integral of an analytic function /(^) is 



fiO = / (/o + M) = h (91) 



and 



j d£,A{i, i) = j di (ao + ai^ + ai^ + 012^^) = ai - 012^ 
j d^A{^, = j di (ao + ai^ + aif + 012^^ = «i + "12^ 
j did^AilO = - j d^diA{lO=-a^2 



(92) 



11 



It is straightforward to show that with these definitions, the following 
expression is a consistent definition of a delta-function: 

5{(,',i)= j dne-'l^^-O (93) 
where ^, ^' and rj are Grassmann variables. 

Finally, given that we have a vector space of analytic functions we can 
define an inner product as follows: 

= j d^d^ e-^^ m 9i0 = /offo + /i<?i (94) 

as expected. 

8.2.3 Properties of Fermion Coherent States 

We defined above the fermion bra and ket coherent states 

m) = e ^ |0), {{^^}\ = {0\e^ (95) 

After a little algebra, using the rules defined above, it is easy to see that the 
following identities hold: 

cMj}) = ^dtel) (96) 

elite-}) = (97) 

(telle. = %(te}l (98) 

(tell 4 = ^.(tell (99) 
The inner product of two coherent states |{^j}) and |{^j|) is 

(tellteD^e ^' (100) 
Similarly, we also have the Resolution of the Identity (which is easy to prove) 

N 

I = j [nt^d^id^i) e i=i |tel)(tell (101) 

Let IV') be some state. Then, we can use Eq. (101) to expand the state 1^/;) in 
fermion coherent states |^), 

N 

IV) = j (nilid^^^Ci) e -1 V(Oltel) (102) 
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where 

^{0 = M^u...,^n) (103) 
We can use the rules derived above to compute the following matrix elements 

(e|c,#) = %^(C) (104) 

{^\c]W = (105) 

which is consistent with what we concluded above. 

Let |0) be the "empty state" (we will not call it the "vacuum" since it is 

not in the sector of the ground state). Let A ^{c]}, {cj}^ be a normal ordered 

operator (with respect to the state |0)). By using the formalism worked out 
above one can show without difficulty that its matrix elements in the coherent 
states 1^) and |^') are 

m ({c]}, {c,}) \0 =ei \ ({C,}, {C^}) (106) 
Thus, for example, the expectation value of the fermion number operator N, 

N = Y,c]ci (107) 



in the coherent state |^) is 

{i\m 



8.2.4 Grassmann Gaussian Integrals 

Let us consider a Gaussian integral over Grassmann variables of the form 

-^[C,C]= j \ \{d^^d(;^ \ e (109) 

where {C,i\ and {C,i] arc a set of 2N Grassmann variables, and the matrix My 
is a complex Hermitian matrix. We will now show that 

Z[C,C] - (detM) e (110) 

Before showing that Eq. (110) is correct let us make a few observations: 

1. Eq. (110) looks like the familiar expression for Gaussian integrals for 
bosons except that instead of a factor of (detM)""^^^ = Pf(M) (here 
PF(M) denotes the pfaffian of M) we get a factor of det M. This is the 
main effect of the statistics!. 
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2. If we had considered a system of N Grassmann variables (instead of 2N) 
we would have obtained instead a factor of VdetM where M would now 
he an N X N real anti-symmetric matrix. 

To prove that Eq. (110) is correct we will consider only the case d = (i — 0, 
since the contribution from these sources is identical to the bosonic case. Using 
the Grassmann identities we can write the exponential factor as 

e =l[{l-^iMij^j) (111) 

The integral that we need to do is 

Z[0, 0]= J (jl d^.^}j n (1 - iiMij^j) (112) 

From the integration rules, wc can easily see that the only non- vanishing terms 
in this expression are those that have the just one and one (for each i). 
Hence we can write 

Z[0, 0] = (- 1)^ j dc}j ^iMi266M236 • • • + permutations 

= (- l)^Mi2M23M34 ■■■ J (j[d^r d^^j • • • ^Jv^jv + permutations 

= (-l)^^Mi2M23M34... Mjv-i, AT + permutations (113) 

What is the contribution of the terms labeled "permutations" ? It is easy to see 
that if we permute any pair of labels, say 2 and 3, we will get a contribution of 
the form 

(-l)2^'^(-l)Afi3M32M24 . . . (114) 

Hence we conclude that the Gaussian Grassmann integral is just the determinant 
of the matrix M, 

Z[0,0]= J n[[d^,dje =detM (115) 

(Alternatively, we can diagonalize the quadratic form and notice that the Jaco- 
bian is "upside-down"). 

8.2.5 Grassmann Path Integrals for Fermions 

We are now ready to give a prescription for the construction of a fermion path 
integral in a general system. Let H the a normal-ordered Hamiltonian (with 
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respect to some reference state |0)) of a system oi fermions. Let be the kct 
at the initial time ti and 1^1//) be the final state at time tf. The matrix element 
of the evolution operator can be written as a Grassmann Path Integral 

-THitf - ti) 

= i VipVij) en X projection operators (116) 

where we have not written down the explicit form of the projection operators 
onto the initial and final states. The action S'(V', V') is 

S{i^, i,) = J^' dt [ihi^dti^ - H{i;, V-)] (117) 

This expression of the fermion path integral holds for any theory of fermions, 
relativistic or not. Notice that it has the same form as the bosonic path integral. 
The only change is that for fermions the determinant appears in the numerator 
while for bosons it is in the denominator! 



8.3 Path integral quantization of Dirac fermions. 

We will now apply the methods we just developed to the case of the Dirac 
Theory. 

Let us define the Dirac field ?Aq(x), with q = 1, . . . , 4. It satisfies the Dirac 
Equation as an equation of motion, 

(i^ - to) V = (118) 

where is a. 4-spinor and = 7'' 9^. Recall that the Dirac 7-matrices satisfy 
the algebra 

{7^7''} = 25'''' (119) 

where g'^'^ is the Minkowski space metric tensor (in the Bjorken-Drell form). 

We saw before that in the quantum field theory description of the Dirac 
theory, ip is an operator acting on the Fock space of (fermionie) states. We also 
saw that the Dirac equation can be regarded as the classical equation of motion 
of the Lagrangian density 

C = 4, {i0 - m) ip (120) 

where ■0 = V'W^- We also noted that the momentum canonically conjugate to 
the field ij} is iV^i from where the standard fermionie equal time anti-commutation 
relations follow 

|VjQ(.x,a;o), V)j(27, 3:o)| = ^ap^^ix - y) (121) 
The Lagrangian density L for a Dirac fermion coupled to sources rja and fia is 

jC = ip {i0 - m) tp + ipT] + fji} (122) 
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and the path-integral is found to be 



= J^i^\Te J |0) 



Vi^Vipe (123) 



where S = J d'^xC. 

From this result it follows that the Dirac propagator is 

iSa0{x-y) = {0\TMx)Mym 



r}=r}=0 



= {x,a\^^^\y,P) (124) 
10 — m 



8.4 Functional determinants. 

We will now discuss more generally how to compute functional determinants. 
Wc have disciisscd before how to do that for path-integrals with a few degrees 
of freedom (i. e. in Quantum Mechanics). We will now generalize these ideas 
to Quantum Field Theory. We will begin by discussing some simple determi- 
nants that show up in systems of fermions and bosons at finite temperature and 
density. 



8.4.1 Functional determinants for Coherent States 

Consider a system of fermions (or bosons) with one-body Hamiltonian h at 
non-zero temperature T and chemical potential /x. The partition function 

Z = tie y J (125) 

H = j dx i){x)'^ h i}{x) (126) 



where /3= l/ksT, 



N ^ dx tpix)'' ip{x) (127) 



and 



is the number operator. Here x denotes both spacial and internal (spin) labels. 
The functional (or path) integral expression for the partition function is 

^ dr tp* (ihdr + h + fj^ ip 

ViP*Vil)e^J ^ ' (128) 

In imaginary time we set t — > —it, with < r < fiti. 



16 



The fields iP{t) can represent cither be bosons, in which case they arc just 
complex functions of x and r, or fermions, in which case they are complex 
Grassmann functions of x and r. The only subtlety resides in the choice of 
boundary conditions 

1. Bosons: 

Since the partition function is a trace, in this case the fields (be complex 
or real) must obey the usual periodic boundary conditions in imaginary 
time, i. e., 

V'(r) = V(r + Ph) (129) 

2. Fermions: 

In the case of fermions the fields are complex Grassmann variables. How- 
ever, if we want to compute a trace it turns out that, due to the anti- 
commutation rules, it is necessary to require the fields to obey anti-periodic 
boundary conditions, i. e. 

t/>(r) = -^(r + Ph) (130) 

Let {|A)} be a complete set of eigenstates of the one-body Hamiltonian h, {e\} 
be its eigenvalue spectrum with A a spectral parameter {i. e. a suitable set of 
quantum numbers spanning the spectrum of h), and {(j^xir)} be the associated 
complete set of eigenfunctions. We now expand the field configurations in the 
basis of eigenfunctions of h, 

i^{T) = Y,i^xMr) (131) 

A 

The eigenfunctions of h are complete and orthonormal. 

Thus if we expand the fields, the path-integral of Eq. (128) becomes (with 
h=l) 

,. / \ - dT^i)l {-dr - £A + m) i>\ 

Z = j m#^dVAj e A (132) 

which becomes (after absorbing all uninteresting constant factors in the inte- 
gration measure) 

Z = Jlpet{-dr-ex+li)Y (133) 

A 

where a — +1 is the results for fermions and a = —1 for bosons. 
Let fpni''') be the solution of eigenvector equation 

(-5, - £A + m) (r) = a„^^ (r) (134) 

where q:„ is the (generally complex) eigenvalue. The eigenfunctions ip^i''') will 
be required to satisfy either periodic or anti-periodic boundary conditions, 

^^{r) = -a^^ir + P) (135) 
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where, once again, cr = ±1. 

The eigenvalue condition, Eq. (134) is solved by 

Mr) = i^n e^""^ (136) 

provided a„ satisfies 

a„ = -iojn + II- ex (137) 
where the Matsubara frequencies are given by (with /cg = 1) 

^^^hnT{n+^), forfermions ^^^^^ 
I 2-jTTn, for bosons 

Let us consider now the function (PaiT) which is an eigenfunction of —dr — sx+n, 

{-dr - ex + fi)fa{T) = a 'fair) (139) 

which satisfies only an initial condition for </5^(0), such as 

^^{0) = 1 (140) 

Notice that since the operator is linear in dr we cannot impose additional con- 
ditions on the derivative of (pa- 
The solution of 

drlmp^ir) = n - e - a (141) 

is 

^o(r) = ^^(0)e('^-^^-")^ (142) 
After imposing the initial condition of Eq. (140), we find 

^A(^) =e-(" + ^^-/^)^ (143) 

But, although this function '/'^(r) satisfies all the requirements, it does not have 
the same zeros as the determinant Det{—dr+fJ. — s\ — a). However, the function 

^F^{T) = l + a^^{r) (144) 

does satisfies all the properties. Indeed, 

„F^{(}) = 1 + a e- + - ^ (145) 

which vanishes for a = a„. Then, a version of Coleman's argument tells us that 

Del {-Or + II - cx - a) 



constant (146) 



where the right hand side is a constant in the sense that it does not depend on 
the choice of the eigenvalues {ex}. 
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Hence, 

Dct {-dr +ii-ex)= const. aFo{/3) (147) 
The partition function is 

Z = e-f^^ = \[ [Det {-dr + M - ^\)Y (148) 

A 

where F is the free energy, which we find it is given by 

F = -aT^lnDet{-dr+ii-ex-a) 

A 

= -aT ^ In (l + ae^(t^ ~ ^^)) + f{Pn) (149) 

A 

which is the correct result for non-interacting fermions and bosons. Here, we 
have set 

f fermions 



-2TA/' In (l - e^/^) bosons 



where Af is the number of states in the spectrum {A}. 

In some cases the spectrum has t he symme try e\ = —S-\, e. g. the Dirac 
theory whose spectrum is e± = i-^/p^ + m?, and these expressions cam be 
simphfied further, 



J\T>ei{-dt + iex) = J{[T)ei{-dt + iex)T>et{-dt- iex)] (151) 

A A>0 

A>0 

after a Wick rotation — > J| Det (-9^ + £^) (152) 

A>0 

This last expression we have encountered before. The result is 

JJ Det {-dl + el) = const. Vo(/3) (153) 

A>0 

where iI}o{t) is the solution of the differential equation 

{-dl + el)Mr)=Q (154) 
which satisfies the initial conditions 

t/'o(0) = drM^) = 1 (155) 

The solution is 

siiih(|rA|r) 

MV = r-n 156 
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Hence 

sinh(|£A|/3) 



MP) = 

|CA| 

.\ex\P 

as /? ^ oo, (157) 



2kA| 

In particular, since 

/^EI-aI -/3E-^ 
]J-^=e ^>o =e ^<o (158) 

we get that the ground state energy Eq is the sum of the single particle energies 
of the occupied (negative energy) states: 



A<0 

8.5 Summciry of the general strategy 

We will now summarize the general strategy to compute path integrals of any 
type. 

We begin with the path integral for the vacuum persistence amplitude (or 
"partition function" ) which we write in general as 

f i [ d'^x £(fields) 

Z V (fields) e J (160) 

and consider a field configuration (pdx) which is a solution of the Classical 
Equation of Motion 

9^1^=0 (161) 
we then expand about this classical solution 

(t){x) = (t)c{,x) + ^{x) (162) 

the action £'[(/)] = J d^x £[</>], 



5[<^e + V] = SicPc] + J d'^xS(t>ix) 



Thus, for example, for the case of a relativistic scalar field with Lagrangian 
density 

jC=l{d^cPf-V{cl>) (164) 
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we write 

£ = £[(/.e] + ^ip{x) [-52 - V"[ct>c]\ ^{x) + ... (165) 
The partition function now takes the form 

Z = Jvcj> e^^I*^] = const. [Det (-9^ - V"[cP,])] e^^I"^-] [1 + . . .] (166) 

For a fermionic system the main difference is that the determinant factor gets 
flipped around 

[Det (diff'erential operator)] — > [Det (differential operator)]"^^ (167) 

where n is the number of real field components. 

Thus, if we wished to compute the ground state energy Eq, we would use 
the fact that 

S[4>a\ = -TEo (168) 

where T is the total time-span and Eq is the ground state energy for the free 
field, and that the partition function for large T becomes 

Z = e~*-^G^, for large real time T (169) 

Hence, 

Eg = lim 1 In^ = ~S[ct>c] + InDet (-9^ - y"[(^e]) + • • • (170) 

1 — »oo 1 1 Z 

Notice that for fermions the sign of the term with the determinant would be the 

opposite. 

In imaginary time we have 

Z = e~(^^G = e-/?-Eo [Ogt (-V^ + ^'[(/.e])] x [1 + . . .] (171) 

and 

Eg = Eo + ^ InDet (-V^ + V"[cj)c]) + . • . (172) 

The correction terms can be computed using perturbation theory; we will do 
this later on. Let us note for now that the way will do this will require the we 
introduce a set of sources ,/ which couple linearly to the field (p, in the form of 
an extra term to the action of the form J dP xJ{x)(j)(x). Upon shifting the fields 
and expanding we get 

iS[(f)c]+i I d^xJ{x)(j)c{x) 
Z[J] = e J 

1/2 'off Ji^)G{x,y)J{y) 
X [Det (-92 - y"[</.e])] ' e^J^Jy (173) 

where G{x, y) is the Green's function, 

G{x,y) = {x\ _^^ _V"[</>e] '^^ ^^^^^ 
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8.6 How to Compute Functional Determinants: 
Heat kernels and (^-function regularization. 

We have seen before that the evaluation of the effects of quantum fluctuations 
involves the calculation of the determinant of a differential operator. In the case 
of non-relativistic single particle Quantum Mechanics wc discussed in detail how 
to calculate a functional determinant of the form Det [—9^ + W{t)] . However 
the method we used for that purpose becomes unmanageably cumbersome if 
applied to the calculation of determinants of partial differential operators of the 
form Det [— + W(x)] , where x = x^. Fortunately there are better and more 
efficient ways of doing such calculations. 

Let A be an operator, and {fn{x)} be a complete set of eigenstates of A, 
with the eigenvalue spectrum S{A) = 

Afn{x) = aMx) (175) 

We will assume that A has a discrete spectrum of real positive eigenvalues. For 
the case of a continuous spectrum wc will put the system in a finite box, which 
makes the spectrum discrete, and take limits at the end of the calculation. 
The function C(s), 

oo ^ 
n— 1 

is the well known Riemann i^-function. We will now use the eigenvalue spectrum 
of the operator A to define the generalized C-fanction 

C4^) = E;^ (177) 



where the sum runs over the labels of the spectrum of the operator A. Then, 
upon differentiation we find 



y 

ds ^ ds 

n 

= -E^ (178) 



Then, in the limit s — > 0+ we find, 



lim ^ = - Vlna„ = -lnTTa„ = -InDetA (179) 
s^o+ ds ^ 



Hence, 



dU 
ds 



= - In Det A (180) 
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Let us define now the generalized heat kernel, 

Ga{x, y; t) = ^ e-«"^ Mx)f:{y) (181) 

n 

where r > 0. The heat kernel Ga{x, y; r) clearly obeys 

-drGA{x, y; r) = AGa{x, y; r) (182) 

which can be regarded as a Heat equation. Indeed, for A = —D V^, this is 
the regular Heat Equation (where D is the diffusion constant); in this case r 
represents time. In general we will refer to r as proper time. 
The heat kernel 0^(0;, ?/;t) satisfies the initial condition 

lim GA{x,y;T) = V fn{x)f*{y) = 6{x - y) (183) 

r— »0+ — 
n 

where we have used the completeness relation of the eigenfunctions {Jn{x)}. 
Hence, Ga{x, y; r) is the solution of a generalized Heat Equation with kernel A. 
It defines a generalized random walk or Markov process. 

We will now show that Ga {x, y; t) is related to the function C,a (s) • Indeed, let 
us consider the heat kernel GA{x,y;T) at short distances, y ^ x, and compute 
the integral (below D is the dimensionality of space-time) 

I d^x \imGA{x,y;T) = e-""-^ f d^'x U{x)f:{x) 

n 

= ^ e-anT = /'tr e'^A (184) 

where we assumed that the eigenfunctions are normahzed to unity 

d'^x \fn{x)f = 1 (185) 



i. e. normalized inside a box. 

We will now use that, for s > 0, 



f 

Jo 



dT e"""^ = ^ (186) 



10 

where r(s) is the Gamma function: 



POO 



no 

r(s) = / dr t'-' e-^ (187) 
Jo 

Then, 

/ rfrr^-i [ d^x lim GA(x,y;T) =Y (188) 
Jo J V < 



23 



Therefore, we find that the generahzed C- function, Ca{s), can be written ob- 
tained from the generahzed heat kernel GA{x,y;T): 



Ca{s) = t^ f dTT' [ d^x \im GA{x,y;T) 
r(s) Jo J y^^ 



(189) 



This resuh suggests the foUowing strategy for the computation of determinants: 

1. Given the Hermitian operator A, wc solve the Generalized Heat Equation 

AGa = -dr Ga (190) 
subject to the initial condition 

lim GA{x,y;T)=6''{x-y) (191) 

T— >0+ 

2. Next we evaluate find the associated function, Ca{s)^ using the expres- 
sion 

U{s) = T^! dr t'-^ [ d^'x \im GA{x,y;T) (192) 

f [s) Jo J v^^ ^ 

tr e~'^^ 

3. We next take the limit s ^ 0+ to relate the C-function to the determinant: 

lim ^%ifl = -lndeti (193) 

s^o+ ds 

In practice we will have to exercise some care in this step since we will find 
singularities as we take this limit. Most often we will keep the points x 
and y apart by a small but finite distance a, which wc will eventually take 
to zero. Hence, we will need to understand in detail the short distance 
behavior of the heat kernel. 

Furthermore, the propagator 

SA{x,y) = {x\A-^\y) (194) 

can also be related to the heat kernel. Indeed by expending Eq.(194) in the 
eigenstates of A, we find 

S{x, y)=Y: ^^^^^ = E ^^^^^^ (195) 

We can now write the following integral of the heat kernel as 

/ dTGAix,y;T) = V /„(x)/:(y) / dr e"""^ 

= Yl ^^^^^^^^ = S{x,y) (196) 
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Hence, the propagator SA{x^y) is an integral over the heat kernel: 

SA{x,y) = / dTGA{x,y;T) (197) 
Jo 

Equivalently, we can say that since Ga satisfies the Heat Equation 

AGA = -drGA =^ GA{x,y;T) = {x\e-^^\y) (198) 

Thus, 

Sa{x, y) = r dT {x\ e-^^ \y) = {x\ i'^ \y) (199) 
Jo 

and that it is indeed the Green's function of A, 

A,S{x,y) = Six-y) (200) 

It is worth to note that the heat kernel Ga{x. y\ t), as can be seen from Eq. 
(198), is also the density matrix of the bounded Hermitian operator A. As such 
it has an imaginary time (r!) path- integral representation. (To actually insure 
convergence we must also require that the spectrum of A be positive.) In that 
picture we view GA{x,y\T) as the amplitude for the imaginary-time (proper 
time) evolution from the initial state \y) to the final state |a;). In other words, 
we picture Sa{x, y) as the amplitude to go from y to a; in an arbitrary time. 



8.6.1 The determinant of the Euclidean Klein-Gordon Operator 

As an example of the use of the heat kernel method we will use it to compute 
the determinant of the Euclidean Klein-Gordon operator. Thus, we will take 
the hermitian operator A to be 

A = -V^ + (201) 

in D Euclidean space-time climensions. This operator has a bounded positive 
spectrum. Here we will be interested in a system with infinite size L ^ oo, and 
a large volume V = L^. We will follow the steps outlined above. 

1. We begin by constructing the heat kernel G{x,y\T). By definition it is 
the solution of the partial differential equation 

(_V2 + m") G{x, y- r) = -drG{x, y; r) (202) 

satisfying the initial condition 

lim G{x,y;T)= 6° {x - y) (203) 

r— >0+ 

We will find G{x, y; r) by Fourier transforms, 

G{x, y;T)=J ^ G{p, r) e 'P' " V) (204) 
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Wc find that in order for G{x,y;T) to satisfy Eq.(202), its Fourier trans- 
form G{p; t) must satisfy the differential equation 



-drG{p;T) = {p^ + m^) G{p;t) 



(205) 



The sohition of this equation, consistent with the initial condition of 
Eq.(203) is 

G{p;T) = e-(p' + '^^)^ (206) 

We can now easily find G{x,y;T) by simply finding the anti- transform of 
G(p;t): 



G{x,y\T) 



f d^p 
J (27r)^ 



[p"^ + m^) T + ip ■ {x - y) 



9 \x — V\ 
At 



[Attt) 



D/2 



(207) 



Notice that for m — > 0, G{x,y;T) reduces to the usual diffusion kernel 
(with unit diffusion constant.) 



\x-y\ 

e 4r 

lim G(x, y; r) = -^-nr- 

— ^ ' {Anrf/^ 



(208) 



2. Next we construct the C- function 

C-v^+„2(s) = -1- / dTT'-' [ d^'x lim G{x,y;T) (209) 

r(s) Jo J 



We first do the integral 



l^drr^-' J d''xG{x,y;T)= dr t^-^-""'^ e 

(210) 

where R= \x — y\. Upon scaling the variable r = \t, with A = R/2m, we 
find that 



2 



Jo 



drr'-^ G{x,y;T) 



{in) 



R 



\2m 



Kn{mR) (211) 



where Kv{z), 



1 r 



(212) 
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is a modified Bessel function. Its short argument behavior is 



-I +... 



(213) 



As a check, wc notice that for s = 1 the integral of Eq.(210) does reproduce 
the Euchdean Klein-Gordon propagator that we discussed earlier in these 
lectures. 



3. The next step is to take the short distance limit 



lim 

fl->o 



f 

Jo 



G{x, y; t) = lim 



21- 



m 



D-2s 



R^o (27r)^/2 (mi?)#-« 



Kn^{mR) 



(214) 



Notice that wc^ havc^ exchanged the order of the limit and the integration. 
Also, after we took the short distance limit i? —> 0, the expression above 
acquired a factor of F (s — -D/2), which is singular as s — D/2 approaches 
zero (or any negative integer). Thus, a small but finite R smears this 
singularity. 

4. Finally we find the ^-function by doing the (trivial) integration over space 

C(s) = / dr j d'^x]iiaG{x,y;T) 

[s) Jo J 



Vis 



T{s- 



(47r)^/2 F(s) 



-2s 



(215) 



where fi = l/R plays the role of a cutoff mass (or momentum) scale that 
we will need to make some quantities dimensionless. The appearance of 
this quantity is also a consequence of the singularities. 



We will now consider the specific case of D = 4 dimensions. For D 
C-function is 



C(s) = V 



-2s 



167r2 (s- l)(s-2) 



4 the 



(216) 



We can now compute the desired (logarithm of the) determinant for = 4 
dimensions: 



InDet r-V^ + m^l = - lim ^ = 

L J „ ,n-t- rto 



m 



0+ ds IGw^ 



In 



m 



V 



where V = L^. A similar calculation for = 2 yields the result 



InDet [-V^ + m^] 



m 
2^ 



In- 



V 



(217) 



(218) 



where V = L . 
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8.7 Path integral for spin. 

We will now discuss the use of path integral methods to describe a quantum 
mechanical spin. Consider a quantum mechanical system which consists of a 
spin in the spin-5 representation of the group SU{2). The space of states of 
the spin-5 representation is 25' + 1-dimcnsional, and it is spanned by the basis 
{\S,M)} which are the eigenstates of the operators and S3, i.e., 

S^\S,M)= S{S + 1)\S,M) 
S3\S,M)= M\S,M) 

(219) 

with \M\ < S (in integer-spaced intervals). This set of states is complete ad 
it forms a basis of this Hilbert space. The operators Si, S2 and S3 obey the 
SU{2) algebra, 

[Sa,Sb]=ieabcSc (220) 

where a,b,c = 1, 2, 3. 

The simplest physical problem involving spin is the coupling to an external 
magnetic field B through the Zeeman interaction 

-f^Zeeman = fJ, B ■ S (221) 

where fj, is the Zeeman coupling constant ( i.e., the product of the Bohr mag- 
neton and the gyromagnetic factor). 

Let us denote by |0) the highest weight state \S,S). Let us define the spin 
raising and lowering operators S^, 

S^=Si± iS2 (222) 

The highest weight state |0) is annihilated by 5+, 

5+|0) =5+|5,5) = (223) 



Clearly, we also have 



5^10)= S'(S' + 1)|0) 
S3\0)= S\0) 



(224) 



Let us consider now the state |n) , 

\n) = e*^(^o >^ " • |0) (225) 

where n is a three-dimensional unit vector {rP = 1), no is a unit vector pointing 
along the direction of the quantization axis {i.e., the "North Pole" of the unit 
sphere) and 9 is the colatitude, (see Fig. 2) 

n • no = cos 9 (226) 
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Figure 1: 



As we will see the state |n) is a coherent spin state which represents a spin 
polarized along the ft axis. The state \n) can be expanded in the basis \S, M), 

s 

l">= E D^Mk^)\S,M) (227) 

M=-S 

Here D)^g{n) are the representation matrices in the spin- 5* representation. 

It is important to note that there are many rotations that lead to the same 
state \n) from the highest weight |0). For example any rotation along the di- 
rection ft results only in a change in the phase of the state \n). These rotations 
are equivalent to a multiplication on the right by a rotation about the z axis. 
However, in Quantum Mechanics this phase has no physically observable conse- 
quence. Hence we will regard all of these states as being physically equivalent. 
In other terms, the states for equivalence classes (or rays) and we must pick 
one and only one state from each class. These rotations are generated by S3, 
the (only) diagonal generator of SU{2). Hence, the physical states are not in 
one-to-one correspondence with the elements of SU (2) but instead with the el- 
ements of the right cosct SU{2)/U{1), with the U{1) generated by S^. (In the 
case of a more general group we must divide out the Maximal Torus generated 
by all the diagonal generators of the group.) In mathematical language, if we 
consider all the rotations at once , the spin coherent states are said to form a 
Hermitian line bundle. 

A consequence of these observations is that the D matrices do not form a 
group under matrix multiplication. Instead they satisfy 

£>(^)(m)D(^)(n2) = D^^\n3) e^*(^i' ^2, r?3)53 (228) 

where the phase factor is usually called a cocycle. Here $(ni,n2,n3) is the 
(oriented) area of the spherical triangle with vertices at ni,n2,n3. However, 
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since the sphere is a closed surface, which area do we actuaUy mean? "Inside" 
or "outsider"? Thus, the phase factor is ambiguous by an amount determined 
by 47r, the total area of the sphere, 

6*47^^ (229) 

However, since M is either an integer or a half-integer this ambiguity in $ has 
no consequence whatsoever, 

e^^^r^ = 1 (230) 
(we can also regard this result as a requirement that M be quantized). 




Figure 2: 

The states |n) are coherent states which satisfy the following properties (see 
Perelomov's book Coherent States). The overlap of two coherent states |ni) and 
1^2) is 

(ni|n2)=(0|I)(^)(^i)^^'^H^2)|0) 

= (0|D(^)(no)e^^^^i'^2,no)53|o) 



(231) 

The (diagonal) matrix element of the spin operator is 

{n\S\n) = Sn (232) 

Finally, the (over-complete) set of coherent states {|n)} have a resolution of the 
identity of the form 

diJ,{n) \n}{n\ (233) 
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where the integration measure rf/i(n) is 

dn{n) = S{n^ - l)d^n (234) 

Let us now use the coherent states {|n)} to find the path integral for a spin. 
In imaginary time r (and with periodic boundary conditions) the path integral 
is simply the partition function 

Z = tve~f^^ (235) 

where /? = l/T (T is the temperature) and H is the Hamiltonian. As usual 
the path integral form of the partition function is found by splitting up the 
imaginary time interval < t < in Nr steps each of length St such that 
NrSr = (3. Hence we have 

Z= lim trfe"*^^-^)^" (236) 

and insert the resolution of the identity at every intermediate time step, 

^=.j^sr^. (n/ ^M".)) (|n(^(-i)i^"'^"''i-(-i+o) 



- N^JZr^o ( n / '^^("^■) ) ( n MT,)|n(r,+i)) - ST{n{Tj)\H\n{Tj+,))] 

(237) 



However, since 

(n(r,)|g|n(r,+i)) 
{n{Tj)\n{Tj+i)) 

and 



{n{rj)\H\niTj)) = i,SB . n{Tj) (238) 



(n(r,)|n(r,+0) = ^ll^^^t^i^y ^^$(r^(r,), n(r,+0, -0)^ (239) 
we can write the partition function in the form 

Z= lim / Vn e^'^^I^] (240) 

N^^oo,St^O J 

where SE[n] is given by 



-SE[n] = iS^^{n{Tj),n{Tj+i),no) 

/-, , ^/ X x\ 



(241) 
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The first term of the r. h. s. ofEq. 251 contains the expression <I>(n(Tj), n(Tj+i), no) 
which has a simple geometric interpretation: it is the sum of the areas of the AT^ 
contiguous spherical triangles. These triangles have the pole no as a common 
vertex, and their other pairs of vertices trace a spherical polygon with vertices at 
{n(rj)}. In the time continuum limit this spherical polygon becomes the history 
of the spin, which traces a closed oriented curve F = {n(T)} (with < t < (3). 
Let us denote by fi"*" the region of the sphere whose boundary is F and which 
contains the pole no- The complement of this region is f2~ and it contains the 
opposite pole —no- Hence we find that 

lim ^{n{Tj),niTj+i),no)=A[n+]=A'K-A[n'] (242) 

iV^— >oo,(5t— >0 

whore A[ii\ is the area of the region fl. Once again, the ambiguity of the area 
leads to the requirement that S should be an integer or a half-integer. 




Figure 3: 

There is a simple an elegant way to write the area enclosed by F. Let n(r) be 
a history and F be the set of points o the 2-sphere traced by n(r) for < r < /3. 
Let us define n(r, s) (with < s < 1) to be an arbitrary extension of n(r) from 
the curve F to the interior of the upper cap f2+, such that 

n(T, 0) = n(T) 
n(T, 1) = no 
n(T,0) = n(T + /3,0) 

(243) 

Then the area can be written in the compact form 

A[n+] = I ds I dr n{T, s) ■ drn{T, s) x d,n{T, s) = Swz[n] (244) 
Jo Jo 
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In Mathematics this expression for the area is caUcd the (simplcctic) 2-forni, 
and in the Physics hterature is usually called a Wess-Zumino action, Swz, or 
Berry's Phase. 




Total Flux 
$ = AttS 



Figure 4: A hairy ball or monopole 
Thus, in the (formal) time continuum limit, the action Se becomes 

Se = -iS Swz[n] + ^ / dr {drn{T)f + dr jiS B ■ n{T) (245) 
2 Jo Jo 

Notice that we have kept (temporarily) a term of order St, which we will drop 
shortly. 

How do we interpret Eq. 245 ? Since n(T) is constrained to be a point on 
the surface of the unit sphere, i.e., = 1, the action S'_E[n] can be interpreted 
as the action of a particle of mass M = SSt and n(r) is the position vector 
of the particle at (imaginary) time r. Thus, the second term is a (vanishingly 
small) kinetic energy term, and the last term of Eq. (245) is a potential energy 
term. 

What is the meaning of the first term? In Eq. (244) we saw that Swz[n], 
the the so-called Wess-Zumino or Berry phase term in the action, is the area of 
the (positively oriented) region ^[0+] "enclosed" by the "path" n(T). In fact, 

Swz[n] = ds drn-drnxdsn (246) 
Jo Jo 

is the area of the oriented surface f7+ whose boundary is the oriented path 
r = 9ri+ (see Fig. 3). Using Stokes Theorem we can write the the expression 
5'^[n] as the circulation of a vector field A[n], 

(f dn ■ A[n{T)] = [[ dS-VftX Anir)] (247) 
J an J Jq+ 
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provided the "magnetic field" x A is "constant" , namely 

B = VnX A[{t)] = S n(r) 
What is the total flux $ of this magnetic field? 



(248) 



dS ■Vn'x A""] 



sphere 

S I dS-n = ^iTS 



Thus, the total number of flux quanta iV^ piercing the unit sphere is 

= — = 2S = magnetic charge 



(249) 



(250) 



We reach the condition that the magnetic charge is quantized, a result known 
as the Dirac quantization condition. 

Is this result consistent with what we know about charged particles in mag- 
netic fields? In particular, how is this result related to the physics of spin? To 
answer these questions we will go back to real time and write the action 



S[ri] 



dt 



M ( dn\ -r'.^, dfi , 

-j +A[nm-^^-,Snit).B 



(251) 



with the constraint n ^ = 1 and where the limit M ^ is implied. 
The classical hamiltonian associated to the action of Eq. (251) is 



H = 



1 

2M 



n X 



[p- A[n]^ ^ + iiSn ■B = Ho + ijlSu ■ B 



It is easy to check that the vector A, 

A = n X {p- 

satisfies the algebra 

[Aa, Afc] = iheabc (Ac - hSuc) 
where a,b,c= 1, 2, 3, Cabc is the (third rank) Lcvi-Civita tensor, and with 

A • n = n • A = 

the generators of rotations for this system are 

L = A + hSn 

The operators L and A satisfy the (joint) algebra 



[La, Lb] = —ihCabcLs 



La, L^ 







(252) 

(253) 
(254) 

(255) 
(256) 

(257) 



[La,nb] = iheabcflc [La,A-b] = iheabcK 
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Hence 



La A' 



^ [La,H]^{) (258) 



since the operators La satisfy the angular momentum algebra, we can diagonal- 
ize and L3 simultaneously. Let \m,t) be the simultaneous eigenstates of 
and L3, 

P\m,l) = h^l{t+l)\m,l) (259) 
Lz\m,t) = hm\m,£) (260) 

where i? = 1 is the radius of the sphere. The eigenvalues £ are of the form 
£ = S + n, \m\ < £, with n e Z+ U {0} and 25 e Z+ U {0}. Hence each level is 
2£ + 1-fold degenerate, or what is equivalent, 2n + 1 + 25-fold degenerate. Then, 
we get 

A^ = L^-n^fi'S^ = L^-h^S^ (262) 

Since M = S6t — » 0, the lowest energy in the spectrum of Hq are those with 
the smallest value of £, i. e. states with n = and £ = S. The degeneracy of 
this "Landau" level is 25* + 1, and the gap to the next excited states diverges 
as M — > 0. Thus, in the M ^ limit, the lowest energy states have the same 
degeneracy as the spin-S representation. Moreover, the operators and L3 
become the corresponding spin operators, thus, the equivalency found is indeed 
correct. 

Thus, wc have shown that the quantum states of a scalar (non-relativistic) 
particle bound to a magnetic monopole of magnetic charge 2S, obeying the Dirac 
quantization condition, are identical to those of those of a spinning particle! 

We close this section with some observations on the semi-classical motion. 
From the (real time) action (already in the M limit) 

S = - [ dt nS n- B + S [ dt [ dsn-dtnxdsfi (263) 
Jo Jo Jo 

we can derive a Classical Equation of Motion by looking at the stationary con- 
figurations. The variation of the second term in Eq. (263) is 

SS = SS [ dt I dsn-dtnxdsn = S [ dtSn{t) ■ n{t) x dtn{t) (264) 
Jo Jo Jo 

the variation of the first term in Eq. (263) is 

i-T rT 



5 

Hence, 



/ dtnSn{t)-B= I dt 5n{t) ■ fiSB (265) 
Jo Jo 

5S = j dt 5n{t) ■ {-ixSB + Sn{t) x dtn{t)^ (266) 
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which impUes that the classical trajectories must satisfy the equation of motion 

liB = nx dtfi (267) 
If we now use the vector identity 

ny. fix dtfi = [fi ■ dtfi) fi — fi ^dtfi (268) 

and 

n • (9tn = 0, and n ^ = 1 (269) 
we get the classical equation of motion 

dtfi = jJiBxfi (270) 

Therefore, the classical motion is precessional with an angular velocity fipr = 
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